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ABSTRACT 


The aberration present when an axial point is imaged by a centered system 
of spherical optical surfaces may be expressed by any one of a number of power 
series. In this paper the reciprocal of the distance between the intersection of a 
ray with the axis and a fixed point on the axis is expressed as a series in h, where h 
is the perpendicular distance from the ray to the fixed point. 

Formulas are derived by which the constant coefficients of the series expressing 
the aberration in the image space of a single spherical surface may be computed 
if the corresponding coefficients for the object space are known, the fixed point 
of reference being the center of curvature of the surface. These formulas can not 
be used in the case of a plane surface since there is no center of curvature. Ac- 
cordingly, after developing transfer formulas by which the aberration may be 
referred to a new point of reference, a second set of formulas is derived in which 
the point of reference is taken as the vertex of the surface. A final set of formulas 
expresses the longitudinal in terms of the 1eciprocal aberration. 

As a numerical example, the computation of the aberration of the third, 
fifth, and seventh orders of an ordinary photographic ‘“‘landscape”’ lens is given 
in full, Its results compare favorably with those of a trigonometric tracing of 
rays in which seven place tables were used. 

In conclusion the convergence of these series is briefly discussed and some rela- 
tions to diffraction theory are pointed out. 
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I. INTRODUCTION 
1. SPHERICAL ABERRATION 


Rays of monochromatic light proceeding from a luminous point 
and passing through a lens are generally focussed not in a single image 


Via Tf 
4 S9¢ 92657A92 94654094 Sac, 


Figure 1.—Spherical aberration of a single lens and of a doublet 


The paths of three rays when the concave lens is not present are shown in full Jines. The con- 
cave Jens and the paths of the rays if both lenses are present are shown in dotted lines. All parts 
of this pare are drawn in proportion to values obtained by an actual calculation. (m1,2~1.54; 
73,4 = 1.63,. 


point, but at varying distances from the lens according to the zone of 
the lens traversed. In the case of an object point on the axis of a 
centered system of spherical surfaces (the plane being considered a 
sphere of infinite radius) this defect in imagery is called ‘‘spherical 
aberration.’’ Points off the axis are affected by the same aberration, 
but they are also affected by other aberrations, the presence of which 
makes analysis more difficult. The term “spherical aberration” is 
sometimes used more broadly to include these extra-axial aberrations 
and sometimes rather loosely in referring to the axial aberration where 
nonspherical surfaces are involyed. In the present paper ‘‘spherical 
aberration” in the stricter sense of the word is discussed and a method 
is given for computing and expressing it quantitatively. 

As an example of spherical aberration consider the portion of Figure 
1 drawn in full lines. Light from an infinitely distant object point is 
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incident upon the lens 44,N;,, and the three rays shown are refracted 
at the two surfaces so as to intersect the axis at S’o4, S/on, and S’2¢, 
which are the foci for these respective zones. If perpendiculars to 
the axis are drawn at these points, their intersections with lines par- 
allel to the axis through the points of incidence on the first surface 
form the locus of the curve S’,P’, which may be considered as repre- 
senting both in direction and magnitude the spherical aberration in 
the image. [or example, the distance Q’,P’s is the measure of the 
longitudinal spherical aberration of the ray at h,=30. The dotted 
lines illustrate the method of correcting the spherical aberration by 
the insertion of a negative component, M43. .NV;, of suitable design. 


2. SERIES EXPRESSIONS OF LONGITUDINAL SPHERICAL ABERRATION 


The aberration of a lens may be expressed quantitatively by giving 
the location of one or more points on the curve or by an equation for 
the curve. Such an equation is customarily written as an infinite 
series in which the terms beyond a certain order are neglected; thus 
the aberration represented by the curve S’,P’, (fig. 1) may be expressed 
as follows: 


Long. spher.=S’—s’=A’h?+ BYV+C/R+ ... (1) 


where S’—s’ represents the abscissa and h the ordinate, and where 
s’, A’, B’, O’,. . . , are constants. The odd powers of h are absent 
because of axial symmetry. 

The number of terms of the series which it is necessary to retain 
depends upon the conditions of the problem in hand. The straight 
line, S’,Q’,, given by the equation S’ =s’ is a first approximation which 
is exact only for infinitesimal apertures. A value of the constant A’ 
may then be found such that the parabola S’=s’+ A’h? gives a suffi- 
ciently close approximation to the curve up to some aperture depending 
upon the accuracy required. Each of the terms of higher order has 
the property of remaining quite small for small apertures and then 
beginning to increase rapidly at some finite value of kh. Accordingly, 
if the accuracy is not to diminish as larger and larger apertures are 
considered, the fifth order! term and then the seventh order! term 
must be taken into account, and so on. 

In the present example A’ is negative while B’ and C’ are positive, 
causing the upper portion of the curve to turn in the positive direction. 

The same aberration, it is evident, could quite as legitimately be 
represented by one of a number of other curves and its corresponding 


1 The longitudinal spherical aberration which varios as the ith power of the aperture (i= 2, 4,6 . . .) corre- 
sponds to an angular aberration (see equation (3), p. 194) which varies as the (i-+-1)th power of the aperture and. 
also to a path difference which varies as the (i-+2)th power of the aperture. ‘Che corresponding reciprocat 
spherical aberration (see Sec. I, 4) varies also as the ith power of tho aperture, but the Iateral aberration 
varies as the (i4+-1)th power. These differences of order of the expressions which refer in ditferent manners 
to the same aberration give rise to a difficulty in nomenclature. To avoid this the use of the terms primary, 
secondary, tertiary, etc., has been proposed. These designations not only become somewhat awkward for 
the higher-order aberrations treated in this paper, but some confusion might occur because primary and 
secondary already have special meanings in connection with astigmatism and curvature of field. 


In the original pokes of L. Seidel the aberrations are measured by their lateral value so that the terms for 
spherical aberration assume tho orders 3, 5, 7, etc, In view of this and the general acceptance of Seidel’s 
work it has been considered as desirable to designate a particular term of spherical aberration by the order 
which the corresponding lateral geometric aberration assumes. ‘The expression “ the third order longitudinal 
(or reciprocal) aberration” may then bo interpreted as referring to “‘the longitudinal (or reciprocal) 
atercebion which corresponds to the third order lateral aberration,” and sitnilarly for the terms of higher 

der. 

It is to be noted that the constant coefficient A(+)) of reciprocal aberration has the dimensions Z-G+) and 
thus corresponds exactly to the above designation of order (t-+1). 
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equation, for instance by a curve the ordinates of which are equal to 
the heights of incidence on the last surface instead of the first. For 
any particular ray the total aberration would necessarily be the same 
according to either equation but, since the two values of h are not 
directly proportional, different portions of the total aberration would 
be attributed to the different orders. It can be shown that the vari- 
ous orders differ, even if only a single surface is considered, according 
as the ordinate is defined as the height of incidence upon the surface, 
as the height of intersection with a plane tangent at the vertex, as the 
distance from the vertex or from the center of curvature of the surface 
measured along a perpendicular to the ray, as the arc, tangent, or 
sine of the angle of inclination of the ray with the axis, or as one of a 
number of other functions of the aperture. Upon investigation it is 
found that the third order coefficient, A’, is identical in each group in 
which (as for the first three cases mentioned) the ordinate reduces to 
the same infinitesimal for very small apertures. For large apertures, 
however, the third order for a particular ray differs in the several 
equations insofar as h? is slightly different. Thus it may be seen 
that, if the expressions “third order,” “fifth order,” etc., are to have 
exact significance, the exact definition of the aperture must be indi- 
cated. Moreover, direct comparisons between different lenses can be 
made only when the aberrations of both are expressed in exactly the 
same manner. 


3. METHODS OF EVALUATING THE CONSTANTS 


It has been pointed out that successively closer approximations to 
the aberration curve can be obtained by evaluating additional terms 
of the series (equation (1)). The usual methods of evaluating these 
coefficients are of three general types. First, if the lens has actually 
been constructed, a number of narrow beams of light may be singled 
out, as in the Hartmann test, by a diaphragm with several small holes, 
and their positions after traversing the lens may be determined 

hotographically or otherwise. Second, from the specifications of the 
ens system the theoretical position of several rays may be found by 
trigonometric ray tracing. In either case an empirical curve or 
equation is fitted to the discrete points thus obtained. Third, 
trigonometric relations, such as those used in ray tracing, are expanded 
as Taylor’s or similar series in ascending powers of whatever function 
of the aperture is chosen as the parameter. The aberration coeffi- 
cients are thereby expressed as functions of the object distance and of 
the constants of the lens system. The present investigation is of this 
third type, which is usually called the algebraic or analytical method. 

Ray tracing gives directly the total longitudinal aberration of one 
or more individual rays for any given lens. The algebraic method, 
however, is often more useful to the lens designer in that it reveals 
more fully the portion of the final aberration which is contributed by 
each surface of the lens and suggests in what manner the design should 
be altered. Moreover, very often when the design is altered only a 
small pate of the computation need be repeated. The algebraic 
method has been considered objectionable because the third order 
equations are not sufficiently exact while if higher order terms are 
included the computation becomes too laborious. It is believed, 
however, that frequently the additional information is worth the labor 
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involved, and, accordingly, formulas are given for enough orders so 
that sufficient accuracy may be attained in practically all cases. 


4. RECIPROCAL SPHERICAL ABERRATION 


The frequent occurrence of reciprocal distances in optical equations 
has suggested that the spherical aberration formulas might be simp- 
lified in form if the aberration instead of being expressed as a differ- 
ence between the focal lengths of different zones were expressed as 
a difference between the powers of the reciprocals of the focal lengths. 

An idea of the relation of the reciprocal aberration to the longi- 
tudinal may be gained as follows: Take for a moment the back 
image distance V,S’,=8", (fig. 1) as the unit of length. Plot a curve 
with the same ordinates as S’,P’,, but with abscissas equal to the 
reciprocals of those of this curve when measured from the vertex, V4, 
as the origin. Where the abscissa of the old curve is slightly less than 
one that of the new is slightly more than one. Thus it is easily seen 
that the new curve will be very similar to S’,P’, but reversed with 
with respect to $’,Q’,. Since the aberration is measured by the 
departure from the vertical straight line, the reciprocal aberration 
in any case is opposite in sign to the longitudinal aberration. The 
statements regarding equations expressing aberration, which were 
made in the two preceding sections (2 and 3), apply directly to the 
new curve except that other symbols, to be introduced later (see 
equations (2) and (2’)) should replace the ones in equation (1). If 
a different unit of length is used the reciprocal aberration curve 
will merely be changed in scale in the horizontal direction. If, 
however, a different origin be selected, the curve will also be changed 
in shape. 


5. HISTORICAL SKETCH AND BIBLIOGRAPHY 


The theory of third-order aberrations is treated in the majority of 
books on geometrical optics. Among these Conrady ? is unexcelled 
from a practical standpoint, while Von Rohr ® gives a more elegant 
ea cao development and also includes an extensive bibli- 
ography. 

STN th-ondar terms for rays in the axial plane were published by 
Keller * and by Bauer.’ Their expressions are not strictly accurate, 
however, in that h is not exactly defined. Kerber® derived the 
fifth order reciprocal aberration of a single surface defining h as the 
height of incidence. Schupmann’ used a fifth order term which 
he credited to Kerber. Von Rohr® and Ké6nig derived the fifth 
order term by Abbe’s method of invariants with the angle of inclina- 
tion of the ray as the variable. Risco ® extended these equations to 
apply to aspherical surfaces. Smart ?° published terms of the fifth 
order, but unfortunately a number of his equations appear to be in 
error. 


14. E. Conrady, Applied Optics and Optical Design, pt, 1, 518 pp., Oxford Uniy. Press, 1920. 

1M. Von Rohr, editor, Die Bilderzeugung in optischen Instrumenten, J. Springer, Berlin, 1904, Also 
translation by R. Kanthack, His Majesty's Stationery Oice, London, 1920. 
aan? A. Keller, Zur Dioptrik, Entwicklung der Glieder fiinfter Ordnung. 24 pp., C. R. Schurich, Munich, 
4K, L. Bauer, Carl's Reportorium f{. Phys. Tech., vol. 1, pp, 219-241, 1836, 
* A. Kerber, Centr. Z. {. Opt. u. Mech., vol. 7, Pp. 217-218, 1886. 

7L. Schupmann, Die Medial-Fernrohre. B. Q. Teubner, Leipzig, 1899. 

1M, Von Rohr, pp. 238-244. Seo footnote 3. 

*M. Martinez R sco, Estudios Gonerales sobre Aberracién Esferica de Orden Superior, Anales Soc. 
Espanola Fis. y Quim vol. 25 pt. 1, pp. 100-138, 1927. 

1° E, H. Smart, Phil. Mag., vol. 20, pp. 82-91, 1910. 
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Early investigators, notably Petzval, may have been in possession 
of formulas of higher orders than the fifth," but it is not known that 
any of these have been published. It is true, however, that investi- 

ations of aberrations from somewhat different standpoints have 
een extended to fifth and higher orders. Most of these are based 
upon the eikonal function, which depends upon the variation of path 
length through the lens. 


6. NOTATION AND SIGN CONVENTION 


A thoroughly consistent and complete system of notation for optical 
equations is, unfortunately, not available. In the present paper the 
aims have been to follow, in general, one of the systems already in 
use and to use as smali a number of symbols as is compatible with 
clarity and brevity in presentation. In particular may be mentioned 
the use of the Greek characters 8 and ¢ with superscripts (6, 6, 
etc.) as constant coefficients of the various orders of reciprocal spher- 
ical aberration. (See equation (2) and footnote 1, p. 189.) 

The original direction of the light is assumed to be from left to 
right. The refracting and reflecting surfaces are numbered in the 
order in which they are encountered by the light and are indicated 
by numerical subscripts. The well-known idea of regarding reflec- 
tion as a special case of refraction in which the index ratio is —1 has 
been extended to include absolute indices (and indices relative to air), 
the index of any medium being considered as negative in sign if the 
light travels through it in a negative direction.’ Nevertheless, any 
other self-consistent sign convention may be used if preferred. 

A list of the characters with their meanings is given herewith. 
Unprimed letters denote magnitudes in the object space, primed 
letters the corresponding magnitude in the image space. If the two 
are identical, as in the case of 7 or ¢, the character is always written 
unprimed. In the following summary the primed character is 
omitted in most cases, it being assumed that its meaning is easily 
deduced from the definition of the corresponding unprimed character. 


ENGLISH LETTERS 


A, B, C,.. . . =the third order, fifth order, seventh order, . . . ., coefficients 
of longitudinal spherical aberration, positive if the aberra- 
tion tends to make the rim ray cross the axis to the right 
of the paraxial focus. 

c=a subscript meaning ‘‘referred to the center of curvature”’ of 
the appropriate surface. 

d=the distance from the old point of reference to the new one to 
which the transfer is made, positive if to the right. 

h=(in Introduction) linear apercure, subject to a variety of 
definitions. 


4 See, for instance, Trans. Opt. Soc. London, vol. 22, p. 214; 1920-21. Also Berek, Grundlagen der prak- 
tischen Optik, p. 43. Walter de Gruyter Co., Berlin, 1930. 2 
4 This convention is also applicable in path length ius Hons as may be seen from the following simple 
example. Suppose a ray from some point A Js reflected back to A from a mirrorat B. The first distance, 
A to B or AB, being positiyo, the distance B to A or BA, after reflection must be negative, Adding the 
distances together would give zero for the path length. However, if the index be considered negative 
after an odd number of reflections and, as fs universally done, the distance in each medium be multiplied 
by the index of refraction to get the equivalent path length in air the result is: 
Equiv, path length=nAB+n’BA 
=nAB+(—n)(—AB) 
=2nAB 


which {s obviously correct, 
Ap alternative convention in which some of the angles are considered as obtuse after reflection has 
been proposed by T, Smith in Trans, Opt. Soc. London, vol. 27, pp, 312-316, 1925-26, 
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h=the distance from the point of reference or origin of coordi- 
nates to a ray, measured along a perpendicular to the ray 
and positive if upward. 

i, j=generalized integers. 

n=the index of refraction of a medium, usually taken relative 
to sir, to be given a negative sign if the light traverses it 
from right to left. 

r=the radius of curvature of » speherical surface, positive if the 
surface is concave to the right. 

s=the distance from the vertex of a surface vo the paraxial object 
point for that surface, positive if to the right. 

t=s—r=the distance from the center of curvature of a surface to 
the paraxial object point for that surface, positive if to the 
right. 

S=the distance from the verlex of a surface to the intersection 
of a ray in the object pencil with the axis; reduces to s as 
the aperture approaches zero. 

T=S—r=the distance from the center of curvature of a surface 
to the intersection of a ray in the object pencil with the 
axis; reduces to ¢ as the aperture approaches zero. 

v=a subscript meaning “referred to the vertex” of the appro- 
priate surface. 


GREEK LETTERS 


a=the acute angle between a ray and the axis, posivive if the ray lies above 
the axis to the left of the point of intersection. 

B@)=the coefficient of reciprocal spherical aberration of the ith order, in the 
series referred to the center of curvature, usually of the sign opposite 
that of the corresponding longitudinal aberration coefficient. 

¢©) =the coefficient of reciprocal spherical aberration of the ith order in the 
series referred to the vertex of the lens surface, usually of the sign 
opposite thac of the corresponding longitudinal aberration coefficient. 

O='he reciprocal of T' (q. v.). 

0, 0’=the angle of incidence, the angle of emergence; that is, the acute angle 
between a ray and the normal to the surface at the point of incidence, 
positive if the ray lies above the normal t9 the right of the surface. 

«=the constant of transfer to a new point of reference. 
=the ratio of the new paraxial reciprocal object (or image) distance to the 
old paraxial reciprocal image distance. 

«,x’ in Part IV take on the special meaning that the transfer is made from the 
vertex to the center of curvature of a single surface, x applying to the 
object pencil and x’ to the image pencil. 

u=relative index=n/n’. 

’=relative index=n’/n. 

p, o, r=the reciprocals of 7, s, and t. 

=the reciprocal of S. 

=the central angle or the angle between the radius drawn to the point of 
incidence and the axis, positive if the radius lies above the axis to the 
left of the center of curvature. 

V=reciprocal object distance from any point of reference on the axis. 


Il. THE RECIPROCAL ABERRATION OF A SINGLE SPHER- 
ICAL SURFACE REFERRED TO THE CENTER OF CURVA- 
TURE 

1. PRELIMINARY STATEMENT OF PROBLEM 


A set of formulas will first be derived which express the reciprocal 
aberration coefficients of the refracted or image pencil of any spherical 
surface when the corresponding constants of the incident or object 
pencil for that surface are known, the point of reference being chosen 
at the center of curvature of the surface. 
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In Figure 2, VS represents the axis of symmetry of a lens system 
of which one surface, with vertex at V and center of curvature at C, 
is shown in cross section. In a symmetrical pencil of rays incident 
upon this surface, consider the ray QP, incident at the point P. Let 
this ray (extended if necessary) intersect the axis at the point S, re- 
moved from the center of curvature by a distance 7’, positive if to 
the right. Let h, denote the perpendicular distance from the center 
of curvature to the ray, positive if upward, and 7 the radius of curva- 
ture of the surface, positive if the surface is concave toward the right. 

Let B®, B®, B™, etc., be constants, such that if hk. is given, T is 
determined to as close an approximation as may be desired by re- 
taining enough terms of the series: 

P=ZTBMREF BONS + BHS+ Wat (1) 

Here ¢ is the limiting value of T as A, approaches zero; that is, the 
paraxial value of 7, 8“? is the third order term of reciprocal aberra- 


V 
—— | T 


Ficure 2.—Refraction of a ray at a single surface 
tion, BA? is the fifth order, ete. For convenience in writing the 
equation, let 7 © and f=, then 
O=1t+BORZ+BOKAT+BMKS+BOHEt ... (2) 


The angles a, 3, and ¢ are shown in eure 2. They are defined by 
the following equations, apparent from the geometry of the figure 


int SOS 
sin a=m=Ohe 


=th,+BPAS+BOKE+ ... (3) 
sin = "= ph, (4) 
g=atd 


where p is written for i. 
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Let the ray pass from a medium of index n into one of index n’ 
: nr’ “ : 
and write for convenience mae p’8 Denote the various functions 


of the refracted ray, PS’, by priming those symbols which when 
unprimed represent the corresponding functions of the incident ray. 
Then: 


C= + POR 2+ ROR At ... (2’) 
sin a’ =7'h’ +p’ h’3+p/OR S+ ... (3’) 
sin 3’= ph’. (4’) 
and 
a +H =d=at9 (5) 


2. DERIVATION OF THE FORMULAS 
By Snell’s law and equations (4) and (4’) 

sin 3=ph.=p' ph’. (6) 

From this h.=y’h’.; then substituting in equation (3): 

sin a=p' th! +p BOR! 3+ pw BOR St pIBOh! T+ p BOR! P+ 

BML hee (7) 
The sines of a, 3, and 3 have now been expressed as functions of 
the same variable, h’,. Of the possible methods of combining these 
to give sin «’, the unknown angle of equation (5), the following one 


using trigonometrical addition formulas and series expansions has 
riers 
been thought the most convenient ;}* 


iss Tho sas 
cos b=1-5 sin? 3 an sinters 
by equation (6) 
= 1 — 5 pp h’ 2-3 pw pth’ a -75 Tedd (i 
a ae Po ee p/10p10f,/ 30 


Leer 1 1 
cos # =1— gPhh’— geht iG ph’ & = 
sin ()— 0’) =sin ¥ cos 8’—cos 8 sin J’ 
1 2 
= (ul ph’ e— 5 (uw) oh’? 
J(u! nlp = Fe Wl — HIP? 
8 16 


i) , / , u , 
— Fag (u" — w/S) ph! P= 55% (nu! — BM) pith’ 


13 The sense of the index ratio and the direction of the ray are thus mutually defined. If the ray passes 
through tho refracting surface from loft to right, then the index n pertains to the medium on the left and yiee 
versa. ‘Tho sign convention for distances and angles is independent of the direction of propagation of 
Rone neta; The index is considered as negative if the light travels from right to left. (See foot- 
note 12, p, 192. 

44 Tt is to be romombered, of course, if ono is interested only {n the third order or the third and fifth orders, 
that the higher orders may bo ignored entirely. Moreover, it will doubtless be found that the method 
of dorivation of the formulas is more easily followed if only the first ono or two terms of each equation are 
takon into consideration, 
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TABLE 1.—Cos 8 cos 0’ 


hit ns h's La 


1 5 7 
puto! —zgu'tet jpn Hap Hel? 


L 1 5 
tutes tyguitet | tgp uises | +55G nite? 


je! + plot +h ult’ +h wtp 
-4 ps tanto +ypguitone 
-% pi +55 p’ipi0 ? 
256 pl 


cos (#—23’)=cos 3 cos 0’ +sin # sin 3 
=the expansion * in Table 1+ p'p*h’? 


=]— , (un! = 1)?p*h’ 2 — 3 (ul? 1)2pth’ 4 


— 7g WIDE) oH’! 
~ 55g (w'2- 1)" ut + 62+ 5) p°h"! 
= os (u/?—1)? (7/5 + 9p/*+ 9p! +7) ph’,1° 
cos a=1— 5 sinta— 5 Ene ee 
=the expansion ' in Table 2 


TABLE 2.—Expansion of cos « 


[The coefficient of each power of A’. is the sum of all the quantities below it in the same column] 


1 he his A’ hs h’ 10 
1 
1| —ge't? | —a'ep =p trB.) —Fu'"B2 =p!'rBO—pBOBO | —p"7B)—p9BBO 
—fulnpion 
1 
—Fultet -4 138 Sutter pity —Fulieerp) —2 y"0798B9 — 
Bo ic 
Suto 
\ 1 3 5 
—jeu"** — ge BO) — 2 no oH a @.a 
— Sut Sue) 
guinea 


sin a’ =sin [a+ (0—0’)] : 

=sin a cos (8— 3’) +cos a sin (8— 9’) 

=the expansion * in Table 3 

=}, +p Oh S +p’ Oh! & +p Oh! 7 4 ae (8) 
the last term being obtained from equation (3’). 


18 This member of the equation has been written in tabular form for clearness and convenience. In Tables 
land 2 the quantities in a single column are to be considered as inclosed in brackets and multiplied by the 
er ofhatthetop. In Tables 3 and 4 the quantity in each small rectangle in the body of the table is to 
multiplied by the power of A at the top of the column and also Dy, the 8—coefficient at the left end of the 
row. T! GUSH, represented by the whole table {s the sum of all these products. The products may be. 
by row or by column (in this case the Jatter is more convenient) and the common factors 
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Since, to the order of approximation attained, this is true for any 
value of h’., the coefficients of like powers of h’, may be equated and 
are readily reduced to the following form: 


r= p'r+(u'—1)p (9a) 
BO = p98) — Yu" pl(u’ —1)%or+ (p+ u!2r?) —p! (p+ =) (9b) 
BO) = p/>BO — You/3BO (u! —1)p[(u’ —1)p+2p’7] 

— ¥en" ol (u’? — 1)? 98x + (p? — p!277)? — p35 (p? — 77)? (9c) 


BO = pM BO — Yy!*BO? (yu! —1)p— Yy’>B)(u! — 1) p[(u’ —1) p+ 2p’7]- 
— Pp? B© p[hh(u!? —1)*p° — p27 (p? — w!47*) + p97 (p?— 77)] 
— Mon" pl (4? —1)?(u’? + 1) p's + (9? = p!21)?(p? + p27?) 
Fe pigt Meta yates 
= p!>(p*—7?)?(p? + 7°)] (9d) 
BPO = BBO — pBOBO (yu! —1)p— Yu!7B (w’ — 1)p[(u! — 1) p+ 2p'7] 
+ Yp!™B © pl (p? — 3/277) — y’ (p? —377)] 
— Yon!®B [dh (u!? — 1)7p8 — w!7(p? — w!47?) + /57(p? —7°)] 
— KuB© plo(u!? — 1)*(u!? + 1) p° = p!?r(p? = !2r?) (6? + 3/274) 
oye! "7, (pa—ra) (pedioae) | 
—Hesu’ p[(m’?—1)?(5u"* + 6p? + 5) p'r + (0? — p7?)?- 
*(5p'+ Gu p?r? + 5y!4r*) — p’7 (9? — 7?) (Spt + 6 p27? + 57*)] (9e) 
B’ (11) oe pg at a 46 4/10/ B(6)2 + 2880 + 7B) (u’ tat l)p 
—Yp"B (u’ — 1) p[(u’ —1)p + 2u’7] 
+ }opB B® pl (p?—3y/*7*) — nw’ (p? —377)] 
= Hou B lh (u!2 — 1)3 09 = p47 (0? = y/277) +b !87 (0? — 77) 
— Hon’™B© pl (ot + 6u!p?7? — 15p/4r4) — y/3(p' + 6 p27? — 1574)] 
— Ke B® p[}6(u’? — 1)? (wu? + 1) p® — w?7(p? — 27?) (0? + 8/27?) 
+ p’®r(p?— 7°) (p? + 37?)] 
—Kou’B plis(u’?—1) (5u!*+ Ou’? + 5) p" — p/27(p?— 277) 
(pt + Qu! ptr? + 5y!4r4) + p!77(p? = 7*) (pt + 2p%7? + 54] 
— Meson’ ol (u’?—1)?(7n’9 + 9n/4 + Ou’? +7) pr 
+ (p?>— p’?7*)? (7 p> + Ou!" pz? + Ou p?rt + 7/57) 
—p'°(p?— 7°)?(7p° + Op'r? + Op?r* + 77°)] (of) 
If, instead, the quantities in Table 3 be factored the following 
equations (written only to the ninth order) result: 


7’ =p" (7+ p)—p or (1’ +p) =n'(7 +p) (10a) 
BO = p98 + Su! (uw! —1) p (e+ 7) (p— u'r) (108) 
Bi = plBO — Fu (wl — 1) plu! = 1)o+ 2u'2] + gn (ul —1)o(o+7) 

(p—p’r) (Ce? + wu’ +1) p?— pw! (u—1) pr + pr" (10c) 
BO = pTBO —Fy!9B(u! — 1)p— u'S8 (u! —1)ol(n! — p+ Qua] 
= an'8 (y! Vela’ +1) (ul? 1) p§ + 4ur(p? + u’r*)] 


1 
+ gH (a! 1) plot 7) (p— n't) (ul + w+ pl + wl + 1) pt 
=p’ (n=l) prt pw? (u’? + 1) iz? = p38 (uw! — 1) p23 + ptr] (10a) 
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(SM STEP Ia SVAN GT Vp 5H78 (w! —Wol(u’—1)p + 2n’7] 

1 2 1 4 7 
— qe BO (u! — 1) ple? + 8n’7"] — gu"*B © (nu! — L)p- 

; 1 , 

[Hl +1) (n= 1) p+ 427 (0° + w’7?)) — FenB (un! — 1)p- 
“[(u! +1) (4-1) p+ 2 p28 (nu! + pw! + 1)p'r + 4y4p?73 + 6y’5r] 
af gn’ =Vpletz)(p— wD) [5Q + pw’ pt p+ p+ pl + I)p* 
=p! (a! —1)(5p* + 6p + 7 yn? + 6p’ + 5) pr 
+ pl? (Sp + p+ 3p? + pl +5)pt7?— p’8(u! — 1) (5p? + 4y’ + 5) p's 
+ p!* (5p! — p’ + 5)p?r*—5p'®(u’ —1)p7° + 5p’?7") (10e) 


If preferred, (7’+u’r) may be substitued for [(u’—1)p+2y’7] 
throughout these equations. In most numerical applications the com- 
puting of the higher powers of y’ may be avoided by dividing each 


equation by the appropriate power of n’, remembering that u iat. 
Equations (9) then take the following form: 


wat+(t-5)p (11a) 


nr nm nn 


oe 82} G-s EN aren 1) -a (B+ os 
pe RP ta f+ G ma ne (11d) 
1 (5) (5) 1 (3) 1 1 
a2 ILCs o[G-s 
ye (eer 
a a nr n\n? 5) | (11¢) 
1 (7) (7) 1 (3)2 1 1 1 (5) 1 1 1 
iat i aro era Ce 7 eLGn wy et 4] 
1g -) 17r/p lr T 
272 lia-z (PE -it( a -5)+4 z(e “)] 
1 
A @-ay Qua )er (BS) ee 
=) p” aes ; 
aa) atm) Go 


If parle light is incident upon a lens surface, r=B =p = 
a» =0. also n=1 then equations (11) are greatly simplified, as 
follows: 
Z 1 
ao (-3)e (12a) 
pg 4 ii\ 
ne -3(1-3) Pay (126) 
78) 1 2 
at-3(I-aa) me (12¢) 


ae -5( ~=n) 6, (12d) 
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The higher-order equations in these last two sets may be written 
in a similar manner. 
3. DISCUSSION 


(a) INHERENT AND PROPAGATED ABERRATION 


If the incident pencil is free from aberration, as is the case when 
it arises from a point on a material object, then B% =B= ... .=0 
and all terms in the aberration of the refracted pencil vanish ex- 
cept those in the upper row of Table 3; that is, all except the last 
term in each of the equations of (9), (10), or (11). This remaining 
term in each order will be referred to as the inherent aberration of 
that order, since it is independent of the aberration due to the other 
surfaces of the system and is entirely due to the surface in question. 
It is, of course, affected by a change of object distance. 

If aberration is present in the incident pencil, then the terms of 
the image aberration which have a @ coefficient do not, in general, 
vanish. These terms are to be considered as representing the aber- 
ration due to the preceding surfaces of the system after its propaga- 
tion through the surface under consideration, and will be desig- 
nated in the following manner: 

The term in the fifth order equation ((9c), (10c), or (11c)) which 
contains the factor 8@ will be referred to as the third order term of 
the fifth order or the p’® (6®) term. Similarly, the third order 
squared term of the seventh order denotes the p’ (8 ©?) term, and 
the third-fifth order term of the ninth order denotes the 6’® (6 
B®) term, etc. These distinctions will apply directly to the system 
of equations to be developed later (Pt. IL) except that the symbol ¢ 
will replace f. 

(b) NOTES ON NUMERICAL APPLICATIONS 


There are a number of similarities among different terms of these 
equations which make the computation easier than would appear at 
first sight. The 6’“(6%) term is identical with the 6’ “**)(g9*) 
term, the B’“*+(69t) term, ete., except for factors in »’*. More- 
over, the 6’“(8)), g/G+2)(B@B®), BEN (g@g_gM™), BG4(g6)2) and 


the 6’“*(8%) terms differ only by factors in »”, 5? and in one 


case tr. (i andj may be 3, 5, 7, or 9; 72j.) 

It is, of course, optional which form of the equations is to be used 
in computing. The third order of (10b) is doubtless simpler than 
that of (9b). However, if orders higher than the third are also to be 
computed, the equations (9) have their advantages, notably in the 
marked similarities among the following bracketed quantities: 


[(u’ = 1)?pr + (p?-+ wr?) — w'(p? + 7) 
[(u’? = 1)?p%x + (p? — 7?) — 1/9 (p? — 7°)?] 
[iu Se EEC Peasy Se Ua “| 


{(u!2 = 1)*(a!2+ 1) 08 + (p= 228) 2(—8-+ wr?) — (0? — 72)? + 7?)] 
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Each of these quantities is the sum of three terms. Each term in 
the last bracket (seventh order) contains as factors the corresponding 
terms in the first two brackets, except for (u/*+1) instead of (u’—1)?. 
Also the terms in the third bracket contain factors of the terms in the 
second bracket. All this reduces considerably the labor involved in 


+ 
+ 
+2 
ie 
te 
p G 
-| 
-2 
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Figure 3.—Zero reciprocal spherical aberration contour 
The point at »’=—1, oat pertains only to orders higher than the third. 
computing. Either of these sets may be written in the reduced 
form, as equations (11). 
(c) ROOTS, ANOMALOUS POINTS, AND SPECIFIC EXAMPLES 


It is interesting to note the conditions under which the aberration 
of a surface vanishes; that is, to find the root of the inherent terms 
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of the equations (9) or (10) when set equal to zero. In equations (10) 
the factors common to all orders may be written thus: 


o[eo (eet) (0) 


RECIPROCAL ABERRATION COEFFICIENTS 


ales =) ° +1 es +3 


RECIPROCAL OBJECT DISTANCE T 


Fiaurs 4.—Reciprocal aberration coefficients for surface of unit 
radius, wp’ =1.6 


When p has a finite value not zero the roots of this are as follows: 


zu’ =0 Physically impossible. 
nu’ =1 Same index on both sides; no refraction. 
tlp=—1 
or Object and image coincide at the vertex. 
iF) 
= 1 


os The well-known aplanatic point. 


T= p/n’ 


These roots are shown in Figure 3. If p=0, the surface is a plane 

and there is no center of curvature from which to measure the dis- 

tances. If the object coincides with the center of curvature then 
139103—32——3 


Sa ee eS ON sr ee A — 
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7=c and the series is nonconvergent. (See Part VII, sec. 1.) 
Another set of formulas will be derived which may be used in these 
cases. 

The orders higher than the third have the additional roots given 
by equating to zero the quantities in brackets in the inherent terms of 
equations (10). The only root yet found for these is the isolated point 


ie} 


+ 1 


RECIPROCAL ABERRATION COEFFICIENTS 


RECIPROCAL OBJECT DISTANCE T 


Ficure 5.—Reciprocal aberration coefficients for surface of unit 
radius, p’=1/1.5 


te 41, y= -1) indicated in Figure 3. Since the value p’=—1 


indicates a reflection, it may be noted that in the image formed by a 
herical mirror the reciprocal aberration of all orders higher than the 
ird vanishes when the object is at a distance 27 from the vertex. 
The advisability of computing and tabulating the values of these 
coefficients for different values of »’ and 7 is being considered. Tig- 
ures 4 and 5 show, for the values p’=1.5 and p’=1/1.5, the form of 
the functions which would result. These values would have to be 
multiplied by the power of p corresponding to the order of the aber- 
ration to obtain the actual value of any particular coefficient. 


8 
t 
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Ill. TRANSFER FORMULAS 
1. DERIVATION 


The refracted or image pencil from one surface is, of course, identical 
with the incident or object pencil of the succeeding surface. How- 
ever, when the point of reference is changed, the aberration coeffi- 
cients are also altered. Formulas will now be developed for com- 
puting the new coefficients. 

Let the subscripts i and j refer, respectively, to a surface whose 
aberration is known and to the succeeding surface, whose aberration 
is to be computed. Then, as before, see (2’) and (3’) for the known 


surface, 
Op =r  tB OM PHB Oh A+... 
sin a! =O gh’ p=7' +P ORB + 2... (13) 
and for the succeeding surface, 
O,=7,+B PAP +BOAS+ .. . (14) 
sin ay=Ojhy=T hy tBO s+ 2. (15) 


ee Bee G ms S 
Ci I Ci 
aes Ly 


Fiaure 6.—Diagram showing aperture, h, as measured from different potnis 
on the axis 


Let C;, and C, (fig. 6) be the respective centers of curyature, and 
let the distance from QC; to C, be d, positive if to the right. Then, 
for any ray, as PS, the adding of axial distances gives: 


1 1 


a7,70, +4 (16) 
Similarly, for a paraxial ray 
1 1 Fs 
a Bd or rae a (17) 
then: 
heltdy= pay = 8 (18) 


where « is a new symbol introduced for convenience. 
Manifestly 
sin a’;=sin ay 


oo 
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or, by equations (13) and (15): 
1 
h’ = a7, O,hy 


Then, substituting in turn from equations (16), (14), and (18): 


h’,=(1+d®,)h; 
=(1+dr,)hj+dB,OR,+ eee 
=Khy+ dB On; +dpOnz+ ... 
This last expression for h’; may be substituted in the right-hand 
member of equation (13), and then since the latter is equal to equa- 
tion (15) 
sin ay=tyhyt BOR, + BOR, + ... 
= the series expanded "* in Table 4. 
Since this equation is true for any value of A, the coefficients of like 
powers of h, may then be equated. 
The first order equation is 
Ty= KT", (18a) 


16 See footnote 15, p. 196. 
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The third order, after combining terms, is 


(1—dr’,)B;® =x°p’, 
or, substituting from equation (18) 
By =xip’,® (185) 
The fifth order, likewise, is 
el —dr,')B;© = 8x7dB; B® + 5B", 
Substituting for (1—d7,’) and 6; from equations (18) and (18d): 
B® =x9(B",) + 3xdp’ ,) (18¢c) 


The higher orders may be reduced in a similar manner and the 
transfer formulas summarized as follows: 


1 

“Toa, v8) 
Ty= kts (18a) 
6, = x89", (184) 
By =x9(B', +3xdB’,?) (18¢c) 
By? = 08(B", + Sud’ (B® + 12x8d?6 ,2) (18d) 

Bs = x98’, +10xdp’ 8’, +5xdp’ ,)? +55x°d?p’ ,229’ 6) 
+55xdp' , 4) (18e) 


By = K!2(67 0) +12«d(p’,B’ +B’ ,©B’,) +78x'd?(p’ 228’. 
+B": B',9) +3640 (98 +273x'd'8’«] (18f) 


2. NOTES ON GENERALITY 


These transfer formulas are more general than has been indicated 
By definition, the series ; 


O=r+PBOR+BORM+ ... 


refers exclusively to the center of curvature of a lens surface. It is 
evident, however, that a similar series may be written referring to any 
point on the axis and that the transfer formulas will apply in trans- 
ferring to any other point on the axis if 0, 7, 6, and h, are replaced by 
the corresponding symbols from the new series. Moreover, these 
formulas apply as well to some of the series with different definitions 
of h mentioned in the introduction, the condition for their validity 


being 
Wh=f(a) 


where ¥ corresponds to © above, but may refer to any point on the 
axis and f(a) is any function of a the use of which results in a con- 
vergent series throughout the required range of h for each case in 
question. 
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3. TRANSFER FORMULAS INVOLVING POWERS OF n 


If the reduced form of the aberration equations (equations (11) 
or (12)) is used, then the transfer formulas must also be divided by 
the corresponding powers of n.. They are then as follows: 


mau ey) 

, as 
et (19a) 
oP aay 
a (is Be de) (19¢) 
=“(Gr 2 as og OR + antares) (194) 


Other shal forms of aberration equations are to be given later. 
The various forms of the transfer formulas are collected here for 
convenience in reference. 


Center of curvature to vertex: 


asap (20) 
rn eo ny 
nyoy=N (Nk) ay (20a) 


7) 


NiO = (nx) (206) 


vei ve +3 (ne) dO cae (20¢) 


gts? = (ne) | BE, a +8 (nx ya eke O12 (ned)? al (20d) 


Vertex to center of curvature: 


= 2 
n n—adni' oy @?) 
1 }} , 
iran ot (21a) 
(3) 4 
Re gaanees (6K) re 13) 9 
LS Ee nit (218) 
p,® 
nitiO43% dinit’ >? | (210) 


ae (sy [x fee m+gs K<d(nit!®) (nit id) 
+12(£d) (nit? | (21d) 
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Vertex to vertex: 
1 


es (22) 
1 =a noe 
Nyoy=KNi/ oy (22a) 
njt{P=Kn lf! ® (226) 
d 
ng Om a8] ni ti 43x (nit! | (22¢) 


N= [ mire + sf (n/t) (’o/) 


+12 (« Vent! oy | (224) 


Higher orders are derived from equations (18e) and (18f) in a simi- 
lar manner. 


IV. THE RECIPROCAL ABERRATION OF A SINGLE SURFACE 
REFERRED TO THE VERTEX 


A second set of formulas will now be investigated in which the point 
of reference is chosen at the vertex of the surface. These are to be 
used in the cases where the formulas referred to the center of curva- 
ture are not applicable (sce p. 203), elthough in many cases either set 
may be used. 

1, DERIVATION OF THE FORMULAS 


As has been stated, aberration series of the form (2) or (2’) may 
be written referring to any point on the axis. Let the two which 
refer to the vertex of a surface (as Vj, fig. 6) be 


DHotf PhP +eEOAs +E MAS + . .. (23) 
for the incident pencil and 
Dao tO Pre On re OR os . ., (237) 


for the refracted pencil. Probably the easiest way to establish a 
relation between these two series is to transfer the aberration to the 
center of curvature by means of equations (18) so that the resulting 
coefficients are identical with those of series (2) and (2’) and then to 
substitute these coefficients in equations (10). 

Let x and x’ be the constants of transfer in the object and image 
spaces, respectively. Then substituting the appropriate values for 
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the object space in the formulas (18), the following relations are 
derived: 


1 
po! trr= PP Et (24) 
T=Ko (24a) 
BP =KE (240) 
BO =8(EO + 3xre%?) (24c) 
Bo = BEM +8xrt ()¢(5) ae 12x?r?¢ @)3) (24d) 


BO = K(C 4 10 «re OE) 4 Ferg O? 455x774 55<arCO) (24e) 


By a similar substitution and also by reason of equations (10a) 
and (24) 


, 
pay MT Save AS AO) Sy 
=— = = = 25 
tn p30 OMe p ior (25) 
in the case of the refracted pencil. Continuing the substitution 
7 =K' 0 =p’ ko’ (25a) 
Bp’ (3) — pls? (3) (25) 


and similarly the higher orders differ in form from (24¢ to 24e) only 
by the use of y’x instead of x. 

Substituting these values throughout equations (10), dividing each 
equation by the power of u’« appearing as a factor of the left member, 


and writing a-3 = 
= (« +2) 
kK 
yO) = yt +2 y(1 yw) 4 eta) ue—o (26d) 
a AK K 


(¢’ (5) +3p/xre! (3)2) 
= n(¢® +84rf ©) — Fut (wb 


[d— 2) 2 +20) +gu(l =») 5(2 +0)(ue— o): 


B 
= 22. 9 
Ke (26a) 


{0 tutu) we +o! (26e) 


(&’ q@) +Sy'xrt’ Ber (5) + 12p/?x?7?¢/ @)3) F 
= we? +-Sxrf SF 6) + 12:7r7¢ 8) —5ug "(1 — up 


ule +80rg*)(1=n)2 (1-4) +20] 
1 : 3 2 
—sut(1— a) 21 tutu? +a) ot duce +40") 


1 4 
+7gu(l = wh(2 +o)(ub = aie tate? +H +n) a 


» 2 p i 
-(i-#') ae +(1 +H) ao (l-n) Ee +o] (26¢) 


bo ty bow vt 


7Y 
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IgA) + 10p’xr¢’ Ge! @) 5p’ «rt! (5)2 


+ bp? rte S)2¢/ (5) + 55y8Ar3¢/ G4) 
=p(t™ + 10crp OEM + Sarg 2 + 55xPr?¢ 8)2¢ (6) 


$+ 55x O#) — (FEO + BerF 2) (1 — yp) —F we 
(EM 4 SxrE OES + 12 azeny(1—n)2| =p) +20 | 
2 1 2 
= Furr —n)o| we +30" |— gus + 3err)(1—y)2- 
ol] (LSP te 5 tau Qotde |= area—p)e- 
bays K? 16 kK 
p? i p! p? 
(lt ppt pt) + 2(ur wit w*) ot 4 a o° + 60° 
1 ‘igeeRER IE : 
+ ypgnt—wS(E +e )( uke) 5 tat at bat at a WS 
5 
— (5+ yt pp —pt—p8)2 ot (St ut Su? + pot 5p!) 5 0? 


3 2 
- (wtp? 5y) So + (5—pt+ 5u°) 5 of — (5 —5z) Pot +508 | 
(26e) 
Equation (266) is used in eliminating ¢’“ from equations (26c) to 


(26e). Following this step, the resulting form of equation (26c) is 
used in eliminating ¢’“ from the higher order equations, and so on. 


Remembering at the same time that, from equation (24)> - =(p—c) 


and ( + °)=. The equations may be brought into the following 


o’=p—p (p—c)=po—(u-1) p (27a) 


1 = pe —F w(u—1)(p—o)*{u(p— 0) —o] (276) 


thd) 7 £45 ul (u-1)(p~o)[(5u+1)(p~ 2) —4e]— 5 u(u-1)- 
-(p—a)? [u(p—o) —o] [(7n?—5u +1) (po)? ~ (54-5) (27¢) 
p-oc)o+o 

£7 = ye —F yE*(u—1)[(Tu + 4)(p— 2) — Ba] 

+5 1 (4-1) (po (74+ 1)(0— 2) ~ bo] 


$5 ne (u— 1)(o—o)[ (633 —33?-7u+1)(p—o)? 
— (96u?—68n—8)(p—«)*a + (40n —32)(p— 0) 0? — 40°] 
16 #4 — DO— 9)*{ulo— 0) — o(83][ut — 4749 + 25u?— 7+ 1)(p+ —0)* 


Le —80p? + 40u—7)(p—o)8o+ (25u*—40n+17)(p—o)?o? 
Be Diovan p—o)a+ (25u*—40n +17) (p—o)?o (274) 
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Ae EON £8(y—1)[(411 w +81)(p—¢)—3300]— pg F®- 


p-o 
“(w= Iu +5)(p~ 0) ~40] + 5 #5 (4 1)(P— o)L Ou +1) (00) 


—80]— 5 ub *(u—1)[ (1984-424? 80u—38)(p— 0)? 

— (240? — 1164-75) (p—o)?o +(70n—51)(p— a) 0? +4107} 

+5 ut ® (u1)(p— 0){ (99: 61u?—9u +1) (po) — (1604? — 1244 
Oye 

(p— 2a + (70-60) (p— a) a? — 624] + gut (u-1)(p— 0): 

-[(429u°— 5314 + 140,38 +20n?-9n +1) (p—o)® 


— (960u*— 1,3723 + 428u2 +38u—10)(p—a)*o 
+(770pu8 — 1,190.2 + 1,320 +20) (p— c)%o?— (280n?—394y +130)- 


-(p—a)?o? + (50n— 40) (p— 7) a*— 40°] — gnu 1)(p—<c)?- 


-[u(p—o) — olf (715y8—1,525y8 +1,335p4— 665p +215p?— 354 +5). 
-(p—a)®— (1,525 u5—3,6714'+2,709)5—1,489u?+3714—35) (p—a)°o 
+ (1,265y'—3,189p8 +2,873yu?— 1,099u +165) (p—«)*o? 

— 575° — 1,369u? +1,0494— 275) (p— a)*o3 + (215? —371n +165)- 
*(p— a)?o*— (354-35) (p— a) 0° +50") (27e) 


It may, at times, be preferable to multiply these equations by n’ 
and write them in the following reduced form 


n'a’ =no—(n—n')p : (28a) 


n!s"® =n¢® —5n(u—1)(p—o)"lu(e— 0) — a] (280) 


nig = nt +5 nb (u—1)(p—2)[(5u+ 1) (pe) —40] 


—En(u-1) (p= 9)"u(p— 0) — ol (Tu? — 5 G86) 


+1)(p—o)?— (54-5) (p—a) a+ 07 


and similarly for higher orders. In computations by means of these 
equations it is not necessary to find the ¢ coeflicients explicitly, but 
only the ‘‘reduced coefficients,” n¢®, n¢™, ete. 


2. DISCUSSION 


It is evident that these formulas are considerably more complex 
in form than those which refer to the center of curvature. They are 
included here as supplementary formulas to be applied where the 
other system is unsuitable. The eleventh order has not been included 
partly because of its unwieldiness and partly because in the cases 
where these equations will be most frequently applied, namely, sur- 
faces of small or zero curvature or where the object is near the center 
of curvature, the inherent aberration will be small and fewer orders 
will, in general, be required than for the other surfaces of the system. 
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The zero points common to all orders of inherent aberration are 


given by 
nu n(E~ 1) | o(2 = 1)+£|e=0 


n=0, physically impossible. 
(u—1)=0, the same root as in the other system. 


2 
(¢- 1) =0,a double root at the center of curvature. 


(u+ (2 = 1)= —1, the aplanatic point again. 


~2 =i (@) +} +2 +3 


Fiaure 7.—Zero reciprocal spherical aberration contour 


The solid lines indicate zero values of all orders, The dotted curyo indicates an 
additional root of tho fifth order. Additional roots of higher orders are not shown. 


These roots are shown for p=1 in Figure 7, as is also an additional 


root of the fifth order. The additional roots of the higher orders have 
not been investigated. 
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Figures 8 and 9, corresponding to Figures 4 and 5, show the values 
of the inherent aberration coefficients for »=1/1.5 and »=1.5 when 
p=1. Itis to be remembered that »=1/1.5 corresponds to yp’ =1.5 in 
Figure 4. 


+4 


+3 


+2) 


wed] 


RECIPROCAL ABERRATION COEFFICIENTS 


RECIPROCAL OBJECT DISTANCE o 
Fiaurs 8.—Reciprocal aberration coefficierts for single surface 
of unit radius, h=1/1.5 


V. INVERSION. TO LONGITUDINAL SPHERICAL 
ABERRATION 


Given the reciprocal aberration of a pencil of rays, the longitudinal 
aberration may be derived by expanding the brackets in the following 
equation 


ah te SIG IDE 
StF TOME... 
56 ite Lee nl 
“a(tt tent 2rOm+ ... ) 
Leona, + } ree 


=s[1+ (sf A2, +f ORM, + sf OH + 8 ORS, + sf H)]-} 


according to the formula (1+2z)"'=1—z2+e¢—-a+ ... If the 
longitudinal aberration be expressed as before. (See Introduction.) 


S=s+ Ah?+ Bh'+ Ch°+ Dh§+ Eh” QQ) 
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the following formulas are found for converting reciprocal into longi- 
tudinal aberration: 


A=—s%® (29a) 

= —$2(¢ ©) — g¢6)2) (296) 
C= =S(FM — 286 EO) + 52¢ G3) (29¢) 
Drea 8 (25 eS 58) ste OE) — gE et) (29d) 


B= —8%{F OV —28(F OE + COEM) + 3g2(F O2EM 4 FE (62) 
—458¢ O3¢ 6) + gi¢ G4] (29e) 


RECIPROCAL ABERRATION COEFFICIENTS 


+2 +3 


RECIPROCAL OBJECT DISTANCE o 


Fiaure 9.—Reciprocal aberration coefficients for single surface of 
unit radius, h=1.5 


These formulas are quite general. If the aberration refers to some 
point other than the vertex of a surface, it is necessary only to replace 
these symbols by the corresponding ones in the new series. The 
parameter may be taken as any function of the aperture and of any 
constants (as, for instance, the constants of a surface of the lens) so 
long as a convergent series result both for the reciprocal and for the 
longitudinal aberrations. 
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VI. NUMERICAL EXAMPLE 
1. DETAILS OF CONSTRUCTION OF THE LENS 


As an illustration of the application of the formulas derived in the 
preceding parts of the paper, the reciprocal spherical aberration of the 
third, fifth, and seventh orders has been computed for an ordinary 
landscape lens such as is found in most low-priced cameras. The 
details of construction of the lens as given by Conrady ” and the 
necessary reciprocals are as follows, except that the decimal point has 
been changed in order to bring it near the first significant figure both 
in the lengths and in their reciprocals: 


r= —0.5556 p= — 1.7999 
M1, 2= 1.63487 di, 2=0.01 nm, 9= 0.61167 

T2= +1.492 po= + 0.67024 
N2, 3= 1.54712 n=, 3= 0.64636 


do, = 0.03 
id Pees 0.2506 P3> — 3.9904 


The formulas (12), (11), and (28) will be used at the respective 
surfaces. The final longitudinal aberration will then be derived and 
compared with the results of trigonometric ray tracing. 


2. FIRST SURFACE; INFINITE OBJECT DISTANCE 


Since 7;=0 and n,;=1, the simplified formulas (12) may be used at 
the first surface. The numerical work is conveniently arranged in 
four columns, one for the paraxial values and one each for the three 
orders of aberration computed, as follows: 


q Seventh 

Paraxial | Third order} Fifth order order 
aa 0.61167 0.61167 0.22885 | 0.08562 
(:- a . 38833 «38833 77115 - 91438 
pis —1.7999 | —1.7999 —5.8310 | —18.890 
wy eres | 1.9816 4. 9816 1, 9816 

4% 6 Mo 

att BO) | 0.60800 | —o.60253 | —1.11381 | —2.1902 


where j = 1, 1, 3, and 5 in the respective columns. 

The first line in this computation need not be written down if the 
computer so prefers. ‘The next four lines contain the factors of the 
coefficients of the respective orders. The reciprocal of the inter- 
section distance measured from the center of curvature could now be 
found for a ray at any height 4’. from the center of curvature after 
refraction at this surface, this quantity being given by the series 


4 3) 7 (5) 4) 
0; =(4 ns ae (Ba \nttht ae a ; Yates + Oe 
1 1 
= —0. 69896 n’,;—0.69253 n’,3h’ 1? — 1.11381 n’15h’a1$ 
—2.1392 n’"h'g®- + + 


al tht a= om 


11 Conrady, Applied Optics, p. $16, See footnote 2, p. 191. 
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3. TRANSFER FROM CENTER OF CURVATURE TO CENTER OF 
CURVATURE; FIRST SURFACE TO SECOND 


Formulas (19), Section III, 3, are used in transferring the point of 
reference to the center of curvature of the second surface. The 
numerical work is most conveniently carried out in two parts, first 
computing the various terms and factors in the formulas, and second 
combining them to find the new coefficients, as follows: 


1 a | +2.0576 x | +0.29840 a | +0.47000 
2 n'a | +3.3639 2 089043 Sats 
3 end | +1.00379 | xt .79237-102 | & uo aU ei aa713s 
4 se . 70600: 10-4 : 
5 (en'd) | +41.00759 | «8 . 62364- 10-4 fo : —. 33214 

Paraxial Third order Fifth order Seventh order 
6 —0. 69253 .11381 | —2. 1392 
! +6. 1912 
8 —4. 0159 
9 +. 0301 

+. 0246-1074 


In this computation the value of d is found from the original data 
(d2=de2—71 +72), then «x and its powers and multiples are com- 

uted, n’d being found incidentally to the computation of x. Line 6 
is copied from the last line of the computation in the preceding sec- 
tion. The ‘“‘sums” in line 9 are the numerical values of the quantities 
in brackets in the formulas, which, when multiplied by the proper 
power of « give the coefficients for the object pencil at the second 
surface. 


4. SECOND SURFACE AND TRANSFER TO THIRD; CENTER OF 
CURVATURE TO VERTEX 


The aberration at the second surface is computed by means of 
formulas (11). As may be seen below, the computation falls readily 
into three parts. In lines 1 to 10 the various monomial and binomial 
quantities and their powers which appear in the formulas are com- 


puted. In lines 11 to 14 are computed the four bracketed quantities 
7 (3) 1 (6) 


which appear in the v —inherent, the nn —inherent, the 


B’ (7) B® p’ (7) 4 : 
Wt (5), and the 7 ~ inherent terms respectively.'® In each 
case the three products composing the quantity are written in the first 
three columns and their sum is written in the last column. In lines 
15 to 20 the various terms in the formulas are evaluated and added 


together to find the constants of the refracted pencil. The inherent 
(3) 


aberration appears in line 16, the terms containing the factor vo are 
in line 17, and so on. 


4 These designations are explained in Sec. II, 3; p. 201. 
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In the transfer to the vertex of the third surface, the first part of 
the computation, lines 25 to 29, has been written below the second 
part, lines 20 to 24 in order to eliminate the recopying of line 20. 

oreoyer, the computation differs slightly: from that in the preceding 
section in the use of formulas (20) instead of (19). 

The numerical work, then, is as follows: 


* 1 1 1 1 
i ps ai Wh dnd 
1 1 +0. 67024 
2 2 +. 44922 
3 3 +. 30109 
4 5 +. 13525 
‘ ‘ ‘ 
(3) | Ge)’ |G) | Gee) +e-0§ 
5 1 —0. 20857 +0, 43322 +0. 40997 —0. 02325 
6 2 +. 043501 +. 18768 +. 16808 —. 23182 
7 SB) [seatatndeaeedse eceecsitecn teas teeaat seine eter —, 44039 
113 ( pin ae 2) 
ne ns nt om non 
8 1 +0, 001203 +0. 14418 +0. 12458 
9 2 +. 001904 +. 020788 ++. 015520 
iets Fee l(a) 
my (St5 (s nm 
10 1 +0, 79192 +0. 23118 +0. 2115S 
Bracketed quantities Sum 
ll B'Q) —0. 00017 +0. 14141 —0. 13676 +0. 00448 
12 B’O) —. 000120 +. 012715 —. 010031 +. 002504 
13 | B’C)(B@)) +. 000143 +. 018394 —. 016795 +. 001742 
4 Bo) —. 0000425 +. 0020395 —. 0021225 +. 0007745 
Paraxial Third order Fifth order Seventh order 
15 = ‘ ae —0. 20857 —0. 54909-10-2 +0. 23329-1073 +0. O246-10-4 
16 | Inher, —. 2148 10-3 —. 3244-10-4 
17 | (8) +. 0231 -10-3 +. O321-10-4 
18 | (B()) —.0119-10-4 
19 | (8(3)2) +. 0035-10-4 
roe +. 0466 -10-1 —. 2761-10~+ 
21 | (x) beac —, 6075 +1074 ++. 1239-104 
22 | (x2) DEE Bers Roar crr eect cto bacccccpoeeoency —. 9275-104 
Ch HET i Wk Se ete bee eect —1, 0797-10-¢ 
24 | noyy nt 1.3520 
Qa 
25 nk a +0. 48888: 10-« 
a Mpg 
27 (ne) 1. 11939103 a =. 08258-10-4 
28 | (n’xd)2 +22, 614 (nx) 1. 18433-10! AEs 
Die —. 75404 | (ne)§ 1, 25303-10¢ | => —. 3418 +1078 
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5. THIRD SURFACE AND INVERSION TO LONGITUDINAL ABERRATION 


The aberration due to the third surface of the lens might also be 
computed by the same method as was used at the second surface. 
However, it is desired to illustrate the use of formulas (28), and more- 
over it is preferable to have the final aberration referred to the back 
vertex without the additional transfer from the center of curvature 
to the vertex. Since the image is in air, then n’ is unity, p=n, 
equations (28) are identical with equations (27), and the ‘‘reduced”’ 
ee n'a’, n’t’®, . . ., are equal to the actual coefficients, 
o,f 

Tn lines 1 to 4 of the computation the preliminary data are pre- 
pared. The bracketed quantities computed in lines 5 to 22 are 
given the designation of the term in which they appear. The vari- 
ous sums of the powers of n appear in the first column of numbers. 
These are multiplied by the appropriate powers of (p—c) and a, 
and the products are written in the next column. Finally, these 
products are added together to give the total value appearing in the 

st column. 

The coefficients of longitudinal spherical aberration are computed 
in lines 28 to 34 by formulas (29), Part V, in much the same manner 
as a transfer is made to a new point of reference. 

These computations are as follows: 


| ] | | | . 
| {| on | (p—0)é of 
| ! |- a “ a= 
| 
| ip 1 | 1.54712 —3, 2304 
2) 2 | 2.39358 +110, 4743 | 
(3) | 3 | 3.70316 —33. 899 
| @ | 4 5.7292 | +109. 711 | 
| | > i 
| Bracketed quantities Sum 
(5) POCO) (6n-+1) 8.7356 -(p—o) | —28, 272 
(6) —4o +3.016 —25, 256 
| 7 | yO(Inh.) (7m. B} 10. 0195 «(p—0)? +104. 947 
(3 —(5n—5 —2. 7356 *(p)o — 6.6759 
(9) | +o? +. 5636 +98, $40 
i—— : = |e 
(10) | CO(EG)2) (7n-+4) 14, 8208 “(p—c) —47. 995 
Qi) | 30 +2, 2621 —45, 733 
| a a Pee es 
ne et A130) (7n-+1) 11.8298 “(p—0) —38, 236 
(13) | —to +4. 524 —33. 762 
G4) | 7) (63n3—. 2 144. 481 “(p—o)3 |—4, 897.8 
15) —(96n2—. —116, 579 AO +020. 75 
} (6) (402-32) 29, 8848 +(p)o? —54.09 
17) } —4o3 +1.71 —4, 030, 3 


t 


v> 
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(is) | ¢/@ (Inh.) (33n'—. .) 65. 0247 *(p—o)* +7, 133. 92 
(19) —(7ni—, % —37. 446 +(p)'or —957. 16 
(20) (25n?—. . 14, 9547 +(p)2o? +89. 06 
(21) —(7n—-7) —3, 8208 +(p) —5.31 
(22) +e! +.32 +6, 260.8 
Paraxial Third order Fifth order | Seventh order 

(23) no; nt; 1. 16660 —0.7827 —0.7709 1.3529 
(24) Inher. +2. 1832 +18. 85+ +4165. 8S +14, 755, 
(25) (¢)) ae e. Sm} —27. 08 —1, 080. 2 
26) (FO) U —35.7 
(27) (F002) -- -- -|---- +119 
(28) Car +1. 0166 +18. 071 +438. 03 +13, 650 
(20) s-term |... _ 21. —15, 573 

« (30) s?-term |. +5,710 
(639) Sum |. .§ +3, 737 
(32, a A By Co ; ++. 98367 —17. 486 —113. 02 —3, 664 
(33) s? 0. 96761 yo yo +7, 915.6 
(34) gr +326. 56 gas +5, 901.3 


Tt is more or less customary to change the dimensions of a lens 
proportionately so that the focal length is 1,000 mm when plotting 
the aberrations. This is done to a satisfactory approximation in the 
present case by merely changing the decimal point. The intersection 
distance measured from the back vertex for an image ray at height 
h’, is then given by the series 


S’5 os 8" + A’h’? + B'h’;! Sr Ch 
= 983.67 — 0.017486 h’,?—0.00011302 h’;4—0.000003664 h’;° 
6. COMPARISON WITH RAY TRACING 


In comparing the results of algebraic and trigonometric computa- 
tions it is well to keep in mind the following essential difference 
between the two: In the trigonometric work the longitudinal spherical 
aberration is found as the difference between the intersection lengths 
of a paraxial and another ray, If either length is determined with 
an uncertainty of 0.01 mm, the uncertainty in the aberration thus 
determined is at least as large as that. In the algebraic series, on the 
other hand, the spherical aberration is expressed as one or more 
correction terms to be applied to the paraxial focal length, and accord- 
ingly the aberration may sometimes be found with an uncertainty of 
0.001 mm or less when the intersection lengths are not known to 
within 0.1 mm. 

For the lens in our example Conrady gives the following values 
obtained by tracing a ray with the use of five place tables: 


hy i 30.0 a= 1.8253° 
s’3= + 983.67 S’,= 967.13 
Long. Spher. = — 16.54 
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For purposes of comparison, the algebraic aberration must be 
determined for the same ray, that is 


h’3=S’3 sin a&’3=967.13 sin 1.8253°=30.805 


This value substituted into the final series in section 5, above, gives: 


A’h’; = —16. 593 s’3= +983. 67 
B’h’;s= —.1018 Long. Spher.= —16. 698 
O'h' f= —.00381 S’,;= +966. 972 


Total = —16. 698 
These results are shown in Figure 10. 


30 


20 


-20. _I5. {for 25: O MM 
LONGITUDINAL SPHERICAL ABERRATION 


Figure 10.—Spherical aberration of “landscape” lens 

The curve shows the aberration of this lens as given by the accompany- 
ing calculation. The circle indicatesa value obtained by trigonomotricray 
tracing using 5placetables. The curve as shown corresponds toa back focal 
length of 983.6 mm. ‘The horizontal scale, however, is twice the vertical. 


In order to determine whether the discrepancy of 0.158 mm 
between the two values of Long. Spher. is due more largely to the 
trigonometric or to the algebraic computation, a seven-figure trigo- 
nometric computation has been carried out. Its results are as follows: 


8’3= 983. 591 a’,=1. 82576° 
S’,= 966. 888 Long. Spher.= — 16. 703 


When this more precise value is taken as a standard of comparison, 
it may be seen that the algebraic computation, even if only the third 
spherical order term be considered, gives a better value for the spheri- 
eal aberration than does the five-figure trigonometric computation. 
Moreover, if the fifth order term is included the resulting aberration 
differs but little from that given by the seven-place work. The sey- 
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enth order term is negligible in the present case, but with larger rela- 
tive apertures it would become more prominent. 

It is interesting to note that Conrady computes the third order 
spherical aberration to be —16.702, which is virtually identical with 
the total aberration given by the 7-place computation, and then 
states that this value “agrees very well with the more exact trigo- 
nometric amount, —16.54, showing that the higher spherical aberra- 
tion is small.’? The difference between the two values of the third 
order aberration may be considered as due to the causes mentioned 
in Part I, section 2. On the other hand, the close agreement between 
the first two paraxial image distances mentioned above appears to 
be largely fortuitous. 

Although a final appraisal of the values of algebraic v. trigonometric 
computations should not be based upon one numerical example alone, 
still the following estimate, based upon this computation along with 
others which the author has made, may well be given here. The 
algebraic method is believed to be capable of giving as dependable 
results as the trigonometric in all practical cases. As a general rule, 
an exactness comparable to that given by 5-figure ray tracing may 
be obtained for apertures up to approximately //10 by computing the 
third order only, and for apertures up to approximately //6 by com- 
puting the third and fifth orders. These stated apertures will vary 
widely, tending to be smaller for meniscus than for symmetrical 
lenses. The labor involved in computing the paraxial and third 
order quantities by the method described in this paper is about the 
same as that of tracing one paraxial and one rim ray, that involved 
in computing also the fifth order is about the same as that of tracing 
one paraxial and two other rays. However, the possibility exists of 
compiling tables of the series aberrations, for instance as functions 
of the index of refraction and the reciprocal object distance, and thus 
of reducing the work to a minimum. 


VII. DISCUSSION 
1. CONVERGENCE OF ABERRATION SERIES 


The validity of series expressions of aberration has been questioned 
on the ground that the magnitude of the terms which are neglected 
is not known. In fact Baker and Filon™ have pointed out two cases 
in which these series actually become nonconvergent. ‘There are 
analogous cases in the system of equations developed in the foregoing 
paragraphs, and these will be briefly discussed here. 

In one case of nonconvergence the image distance for rays at some 
finite distance from the axis is infinite and, accordingly, the longitu- 
dinal aberration is infinite and can be expressed by no convergent 
series. In reciprocal aberration this case can be seen merely to be 
one of a zero value of the series. Nevertheless, analagous difficulty 
arises in the reciprocal series when a ray passes through the point of 
reference so that the reciprocal distance, and hence the reciprocal 
aberration, becomes infinite. This difficulty, however, can easily be 
avoided by a proper choice of the point of reference. 

The second case mentioned by Baker and Filon concerns aberration 
series which are written in terms of the inclination of the ray in the 


1" Baker and Filon, On an Empirical Formula for the Longitudinal Spherical Aberration of a Thick 
Lens; Trans, Opt. Soc. London, vol. £0, pp. 67-92, 1918-19. See particularly pp. 74-78. 
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image space and involves longitudinal and reciprocal series alike. It 
was shown that two rays in the image pencil may have equal inclina- 
tions but unequal aberrations so that the aberration is a double- 
valued function of the independent variable. Manifestly no series 
can be convergent for both values. 

In the present system of equations there is the following analagous 
case: Let ST (fig. 11) be the caustic curve of a pencil of rays. If the 
origin be chosen as at O so that a normal, N, to the caustic curve 
passes through it, then for every ray, RpSp, tangent to the curve 
bevond N, there is another ray, RaS,, such that the normals, OQ and 
OP, of the two rays are equal. There are, then two values of the 
aberration for a single value of the variable, h, and consequently no 
series can be convergent throughout this interval. 


EG O SPSS 
Ficgurr 11.—Diagram illustrating nonconvergence of aberration series 


However, by choosing a suitable point of reference, this difficulty is 
avoided, or at least the interval of convergence may be sufliciently 
extended for practical purposes. 

Without stopping for the proof for other types of series, it may be 
stated as a more or less general rule that, if any aberration series be 
written using as a variable a length measured from a point on the axis 
to a point on the ray in some specified manner, than an absolute limit 
to the interval of convergence of the series is reached when the speci- 
fied point on the ray coincides with the point of tangency of the ray 
with the caustic curve. In practice, however, long before this limit 
is reached the series will become so slowly convergent as to be of no 
practical yalue. It is to be observed, however, that nonconvergence 
can be avoided in the proposed system by a proper choice of the point 
of reference, while Baker and Filon found it necessary to use a func- 
tion of the incident ray rather than of the refracted ray as the variable. 

This is objectionable on account of the indefiniteness inherent in 
using a function of the incident ray. For example, if A is given in 
SP, (fig. 1), S’; may be found from equation (1). However, this 
determines only one point on the ray since the ray might traverse the 
lens by a different path and still intersect the axis at the same point. 
Even if the height of incidence on the last surface is used as the ordi- 
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nate and the axial point determined as before, the exact direction of 
the ray is not known unless the radius of the surface is given also. 
On the other hand, in the proposed system when the axial point is 
found by equation (1), the position of the ray is competely determined 
by a simple geometrical construction. Moreover, the ratio (or in 
the case of reciprocal aberration, the product) of these two quantities 
equals the sine of the angle of inclination of the ray while in the two 
systems just mentioned there is no direct relation. 


2. RELATION TO PATH DIFFERENCE EQUATIONS 


If it is desired to find the path difference between the two rays at 
some point on the axis, in order to work out the diffraction effects 
according to formulas such as those published by L. C. Martin,*? a 
further advantage is gained in measuring the aperture along a line 
perpendicular to the ray rather than to the axis. In Martin’s formula 
the independent variable is the height of intersection of the ray with a 
wave front through a fixed axial point. If HR4R, (fig. 11) is a wave 
front, then R4G is such a height. Plainly an axial aberration series 
could not easily be written using this variable since the shape of the 
wave front is not usually known. If the wave front were exactly 
spherical then it is evident that 224@ would be exactly equal to EF 
which is the variable used in the present paper. Now, since wave 
fronts do not ordinarily depart from sphericity by more than one or 
two wave lengths, or 0.001 mm, when equals, say, 20 mm, the two 
lengths could be considered equal in most cases without appreciable 
error, so that the coeflicients of the series derived in this paper could 
be substituted directly into Martin’s eqnations. 

A further investigation of the relation between reciprocal aberration 
and path length differences has been begun, and exact equations have 
been derived. It is expected that these will be published soon. 


3. CONCLUSION 


Besides the theoretical interest which attaches to the equations 
derived in this paper, there also is believed to be a considerable 
possibility of their practical application. No optical computer to-day 
questions the value of the third order aberrations developed by Seidel, 
and yet Berek” states that ‘“L. Seidel himself did not consider the 
aberration formulas of the third order which he published between 
1853-1856 suitable for practical purposes and, therefore, furnished 
trigonometrical tracing formulas to the C. A. Steinheil Optical Co. 
(1866).”’ Further, it may be stated that the mathematical tools 
furnished by Seidel lay idle for about 60 years, and only recently 
their usefulness, in spite of their restrictions, has been recognized. 
Although the fifth order, up to the present time, has been considered 
too unwieldy to be useful, still it is not impossible that it, like the 
third order, will in time be brought into a simpler form adaptable to 
numerical work. In addition to this there is the possibility of the 
compilation of tables of these aberrations which would shorten 
considerably the work of computing. 


WasHINGTON, June 16, 1932. 


aon Eee Martin, A Physical Study of Spherical Aberration, Trans. Opt. Soc,, London, vol. 23, pp. 63-92, 
au Berek, p. 41, see footnote 11, p. 192. 
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